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Abstract 

Within the framework of three-neutrino and four-neutrino scenarios that 
can describe the results of the LSND experiment, we consider the capabilities 
of short baseline neutrino oscillation experiments at a neutrino factory. We 
find that, when short baseline (L < 100 km) neutrino factory measurements 
are used together with other accelerator-based oscillation results, the complete 
three-neutrino parameter space can best be determined by measuring the rate 
of v e — > v T oscillations, and measuring CP violation with either v e — > i>n 
or Va — > v T oscillations (including the corresponding antineutrino channels). 
With measurements of CP violation in both v e — > and — > v T it may be 
possible to distinguish between the three- and four-neutrino cases. 



I. INTRODUCTION 

The recent results from the Super-Kamiokande (SuperK) detector confirm the solar [l 
and atmospheric neutrino deficits and strongly suggest the existence of neutrino oscil- 
lations. The azimuthal angle and energy dependence of the atmospheric data indicates a 
mass-squared difference scale, £m 2 tm ~ 3.5 x 10~ 3 eV 2 . The LSND measurements || indi- 
cate neutrino oscillations with a different scale, Sm1 SND ~ 0.3-2.0 eV 2 ; there is also a small 
region of acceptable parameter space at 6 eV 2 . The evidence for oscillations in solar neutrino 
data [@],|], when taken as a whole, prefer yet a third, lower mass-squared difference scale, 
<5 m sun ^ 10~ 4 e V 2 - Once oscillation explanations for some or all of the data are accepted, 
the next step is to attempt to find a neutrino mass and mixing pattern that can provide a 
unified description of all the relevant neutrino data. 



In this paper we consider two possible neutrino scenarios which have been proposed to 
account for the LSND results: (i) a three-neutrino model that also describes the atmospheric 
neutrino data (in which case the solar data would be explained by a phenomenon other 
than oscillations 0), and (ii) a four- neutrino model that also describes both the solar and 
atmospheric data ||[7j . As we show, the three-neutrino model may be considered a sub-case 
of the four-neutrino model. We examine the ability of short-baseline experiments at a muon 



storage ring-based neutrino factory [g-uj to determine the oscillation parameters in each 
case, and discuss how the three- and four-neutrino scenarios may be distinguished. 

We concentrate on experiments in which muon and tau neutrinos are detected via charge- 
current interactions and there is good sign determination of the detected muons and tau 
leptons. The sign determination allows one to distinguish between v e — > and — > z/ M 
for stored \i~ in the ring, and between v e — > and — > for stored ji + . We consider 
primarily the appearance channels v e — > u^, v e — > u T , and — > v T (and the corresponding 
antineutrino channels), for which the uncertainty on the beam flux is not critical to the 
sensitivity of the measurements. We find that the — > v T and u e — > channels provide 
the best sensitivity to the CP- violating phase; the v e — > v T channel allows measurement of an 
independent combination of mixing parameters. Short-baseline experiments at a neutrino 
factory are sufficient to determine the complete three-neutrino parameter space in these 
scenarios when their results are combined with other accelerator-based experiments that are 
currently underway or being planned. Measuring CP- violation in both — > v T and v e — + 
allows one to possibly distinguish between the three- and four-neutrino cases. Hence all three 
off-diagonal oscillation channels are useful. In the four-neutrino case it is not possible to 
have complete mixing and 5m 2 parameter determinations without additional measurements 
in future long-baseline experiments. 



II. THREE-NEUTRINO MODELS 

Here we address the ability of short-baseline neutrino oscillation experiments to probe 
a class of three-neutrino models that can describe the results of the LSND experiment, 
together with the atmospheric neutrino deficit observed by the SuperK experiment. 



A. Oscillation Formalism 

In a three-neutrino model the neutrino flavor eigenstates v a are related to the mass 
eigenstates Vj in vacuum by a unitary matrix U, 

\v a ) = J2 U *j > (!) 

3 

with a = e, fi,T and j = 1, 2, 3. The Maki-Nakagawa-Sakata (MNS) flTj mixing matrix can 
be parameterized by 

(C13C12 C13S12 s 13 e lS \ 

-C23S12 - Si3S23Cl2e i<5 C23C12 - S 13 S 2 3S 12 e lS C13S23 , (2) 
iS i5 j 

^23 s 12 — Si3C23Ci2e — S23C12 — 513023^126 C13C23 / 
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where Cjk = cos Sjk = sin and 5 is the CP non- conserving phase. Two additional 
diagonal phases are present in U for Majorana neutrinos, but these do not affect oscillation 
probabilities. 

For three-neutrino models in which the results of the atmospheric and LSND experiments 
are explained by neutrino oscillations there are two independent mass-squared differences: 
^ m LSND = 0.3-2.0 eV 2 and £m 2 tm ~ 3.5 x 1CT 3 eV 2 . We take 5m\ 2 = 5m 2 tm and bm\ x = 
5m1 SND ; see Fig. |I]a. These mass-squared differences obey the condition Sm1 SND ^> Sml tm . 
Then the general off-diagonal vacuum oscillation probability is 

P(u a -> v p ) ~ A\U al \ 2 \Upi\ 2 sin 2 A L snd - ARe(U a2 U* a3 U; 2 U p3 ) sin 2 A atm ± 2 JS , (3) 

where 

Aj = 1 .275m 2 (eV 2 )L(km)/^(GeV) , (4) 
S = [sin2A atm + sin2A LSND - sin2(A LSND + A atm )] 

= 2(sin 2A atm sin 2 A L snd + sin 2A LS nd sin 2 A atm ) , (5) 



and J is the CP-violating invariant [|H|[nj, which can be defined as J = Im{U e 2U* 3 U* 2 U fJ , s }. 
The plus (minus) sign in Eq. (||D is used when a and j3 are in cyclic (anticyclic) order, where 
cyclic order is defined as e/rr. For antineutrinos, the sign of the CP-violating term is 
reversed. For the mixing matrix in Eq. (|^), 

J = Si 3 C 2 13 S 12 C 12 S 23 C 23 SU1 5 . (6) 

For recent discussions of CP violation in neutrino oscillations, see Refs. [PI, [Ol- |bl and 



22| . For L < fOO km, matter effects |23| are very small and the vacuum formulas are 
a good approximation to the true oscillation probabilities. 

The class of scenarios that we are considering is designed to account for the large — > v T 
mixing of atmospheric neutrinos and the small — > v e mixing in the LSND experiment. 
The v e survival probability at the leading oscillation scale is given by 

P(u e -> u e ) = P{v e -> u e ) = 4c 2 2 c 2 3 (l - c\ 2 c\ 3 ) sin 2 A LS nd • (7) 

Results from the BUGEY reactor experiment p4| put an upper bound on P[p e — > v e ) for 
^ m LSND > 0.01 eV 2 , which provides the approximate constraint 

s\ 2 + s 2 13 < 0.01 . (8) 

This leads to the conditions 

012, 013 «C #23 , (9) 

with the mixing of atmospheric neutrinos, #23, near maximal (#23 ~ 7r /4). Hence we can 
take Si2,si3 <C S23,C23 and c\ 2 — C13 ~ 1. The off-diagonal oscillation probabilities are (to 
leading order in the small mixing angle parameters) 

P{y e -> ~ 4|si 2 c 23 + s 13 s 23 e tS \ 2 sin 2 A LSND - 2si2Si 3 c,5 sin 26> 23 sin 2 A atm 

+s 12 s 13 s s sin 29 23 S , (10) 

P(u e -> u T ) ~ 4|si 2 s 23 - s 13 C23e J<5 | 2 sin 2 A LS nd + 2si2Si3C,5 sin 26» 23 sin 2 A atm 

-Si 2 s 13 s & sin26 23 S , (11) 

P(^ -> v T ) — 4|si2C 2 3 + Si 3 s 2 3e ll5 | 2 |si2S23 - Si 3 c 23 e" 4<5 | 2 sin 2 A LSN d 

+ sin 2 26*23 sin 2 A atm + Si 2 s V3 s 5 sin 29 23 S , (12) 
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where eg = cos 5 and ss = sin 5. The expressions for the antineutrino channels are obtained 



by changing the sign of the .S-term in each case. A representative scenario, 1B1 in Ref. p5 
that we will use as an example has the parameters 

5ml 2 = 3.5 x KT 3 eV 2 , sin 2 29 23 = 1.0 

Srr^i = 0.3 eV 2 , sin 2 29 12 = sin 2 29 13 = 0.015 , (13) 

with 5 a varied parameter. More generally, <5^ 2 snd m the range 0.3-2.0 eV 2 is allowed with 
f/u — * v e oscillation amplitude given by 



tfi2 /0.06eV 2 V 
4|si2C23 + si3S2 3 e J T - "T^ • ( 14 ) 

\ "^"LSND / 



In a short-baseline experiment, L/E should optimally be chosen such that Alsnd ~ 1) 
in which case A atm ~ Alsnd/100 <C 1; e.g., for <5m 2 SND = 0.3 eV 2 and E v = 14 GeV we 
might choose L = 45 km, giving: 

A atm = 0.014 ( 5m \ 2 ) (L) ( li^X) . (15) 
^3.5 x 10- 3 eV 2 J V45 km A E v ) K J 

Note that in the forward direction E v = 14 GeV is the average energy for stored unpo- 
larized fir with = 20 GeV. Thus in the probability equations above, su, si 3 , and A atm 
are all small parameters, at the few percent level or less. Then to leading order in A atm 

S ~ 4A atm sin 2 A L snd • (16) 

Therefore, in these scenarios, the dominant contribution to P{y e — ► v^) and P{y e — > v T ) 
comes from the leading oscillation (the sin 2 Alsnd term), but the dominant contribution 
to P(V(j, — > u T ) comes from the subleading oscillation (the sin 2 A atm term); these results 
are summarized in Table Q. In each case, the dominant term is CP-conserving and pro- 
portional to the product of two small parameters. The CP- violating contribution in each 
case is 2s 12 Si 3 S5 sin 26* 2 3 sin 2A atm sin 2 A L snd, which is a product of three small parameters 
(assuming 5 is not small), and hence is smaller than the CP-conserving contribution (see 
Table |). 

TABLE I. Leading and next-to- leading contributions to neutrino oscillations in the 
three- neutrino models under consideration in a short-baseline experiment with L/E chosen such 
that Alsnd ~ 1- The parameters S12, si 3 and A atm ~ Alsnd/100 are all <C 1. Also shown is the 
size of the CP-violating [CPV) term compared to the (dominant) CP-conserving (CPC) term, 
assuming terms the same order of magnitude cancel in the ratio. 



Oscillat 


ion 


Leading term 




Next-to-leading term 


Ratio of CPV 






(CP conserving) 




(CP-violating) 


to CPC terms 


P{v e ^ 




4|si 2 C23 + si3S23e ,,5 | 2 sm 2 


Alsnd 


4si2Si3S,5 sin 26*23 sin 2 A L snd A atm 


S<5A a tm 


P(Ve^ 


v T ) 


4 S12S23 - si3C 2 3e lS 2 sin 2 


Alsnd 


-4si2Si 3 s,5 sin 26*23 sin 2 A L sNDA at m 


-S^Aatm 


P(^^ 




sin 2 20 23 A 2 tm 




^si 2 s 13 s s sin 26*23 sin 2 A L snd A atm 
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In all, there are six parameters to be determined in the expressions in Eqs. (|10|)-(|T2|): 
the three mixing angles, the phase 5, and two independent mass-squared differences. How- 
ever, although there are six off-diagonal measurements possible with \i and r detection 
[the three in Eqs. (|T0|)-(|T2|) plus the corresponding antineutrino channels], the leading and 
next-to-leading terms in the expressions for these oscillation probabilities are only sensitive 
to five independent quantities: \si 2 c 2 z + si 3 S23e Jl5 |, |si 2 s 2 3 - si 3 c 2 3e i<5 |, Si 2 s 13 ss, sin 2 A LS nd, 
and sin26 , 23A atm . For the parameter ranges we are considering, #12, #13 > A atm ; then the 
sin 2 Alsnd term in Eq. fll2|) can be comparable to one of the other terms in Piy^ — > v T ), 
but this term still depends on a subset of these same five independent quantities. 

The K2K J2§, MINOS J27|, ICANOE |§ and OPERA J2T§ long-baseline experiments 



will measure the parameters in the leading term of the z/ M — > v T probability, #23 and <5m 2 
and MiniBooNE |30[ will measure the parameters in the leading term of the v e — > prob 



atm 1 



ability, |si 2 c 2 3 + Si 3 s 2 3e i,5 | and Sm1 SND . Therefore, only two independent quantities will 
remain to be measured in short baseline experiments: (i) the amplitude of the leading os- 
cillation in the v e — > v T probability, and (ii) the subleading CPV term, which has the same 
magnitude for each off-diagonal channel and can be determined by a comparison of neutrino 
to antineutrino rates. Hence a combination of short-baseline measurements with the results 
of the other accelerator-based experiments would allow all of the parameters in this three- 
neutrino scenario to be determined. Measurements of the other short-baseline off-diagonal 
oscillation probabilities may then be used to check the consistency of the result and/or 
improve the accuracy of the parameter determinations. 

One can also measure the — > survival probability, which to leading order in small 
quantities can be written 

P(z/ M Vfl ) ~ 1 - 4A(1 - A) sin 2 A L snd - sin 2 26 23 sin 2 A atm , (17) 

where 

A — S 23 S 13 + C 23 S 12 + 2S23C23Si2Si 3 C,5 . (18) 



In short-baseline experiments, both oscillatory terms in Eq. ( |TTD are second order in small 
quantities. A measurement of P{v^ — > v^) could be used in conjunction with other short- 
baseline measurements to make a complete determination of the oscillation parameters with- 
out the use of other data. However, since A and A atm are small, the deviations of Piy^ — > v^) 
from unity are also small, and the normalization of the beam flux would need to be known 
to high precision for this to be a useful measurement. 

An examination of Eqs. (|T0|)-(|T2D and (|16|) shows that only the relative sign of 5 and 



5m\ 2 may be determined in short-baseline measurements for the distances and oscillation 
parameter values that we are considering. Long-baseline experiments at a neutrino factory 
should be able to determine the sign of 5m\ 2 since there is a significant dependence of the 



matter effect on the sign of 5m^ 2 for L > 2000 km [12,13,15 . 



B. Results 



Figure |2] shows the ratio of the rate of tau production for stored fi + to that for stored 
for the v e —>■ v T and —>■ v T channels, versus the CP phase 5 for several values of 
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baseline length L, with oscillation parameters given by Eq. flT3|). Statistical errors are also 
shown, assuming 10 20 kt-decays (corresponding, for example, to three years of running with 
10 20 useful muon decays per year and a 1 kt detector having 33% tau detection efficiency). 
The event rate calculations are done according to the method outlined in Ref . |ll| • F° r the 
choice of parameters in Eq. (|13"D, s 12 and S13 are larger than A atm and the — > v T channel 
shows the largest relative CP- violating effect (see the last column of Table |). 

The last column in Table [I] shows that the relative size of the CP- violating effect in the 
Vy, — > v T channel decreases with increasing A atm (i.e., with increasing L/E). Also, for very 
small L/E, when A L snd < 1, 5 -> to leading order in A atm and A L snd, and the CP 
violation becomes negligible. Therefore, for any given set of oscillation parameters, there 
will be an optimum L/E that maximizes the CP violation effects in the — > v T channel. 

Figure ||a shows the ratio of v^, — > v T event rates (from /i + decays) to — > v T event rates 
(from yT decays), P MT , versus baseline for 20 GeV muons, for oscillation parameters given 
by Eq. (|13|) and three values of 5 (90°, 0°, and —90°). Approximate statistical errors are 
shown for 10 20 kt-decays. The figure clearly shows that there is one distance that maximizes 
the size of the CP- violation effect, which in this case (20 GeV muons) is about L = 45 km. 
This optimal distance decreases slowly with increasing 5m| 2 . For 8rn\ 2 in the range 2.5- 
4.5 x 10~ 3 eV 2 , we find that the optimal L is in the range 40-50 km for <5m 2 SND = 0.3 eV 2 . 
The optimal L scales inversely with <5mL SND , and for stored muon energies well above the 
tau threshold may be approximated by 

e.g., for <5^lsnd = ^ e ^ 2 ^ ne ^° es ^ seris itivity to CP violation is obtained with L ~ 6 km. The 
distance from Fermilab to Argonne is about 30 km, which would be optimal for E^ = 20 GeV 
and (^lsnd = 0-45 eV 2 . Similar results for an optimal distance for CP violation in — > v T 
in the context of four- neutrino models have been reported in Ref. [f2~T|]. 

Given Eq. ([14]) (the LSND constraint on the — > v e oscillation amplitude), CP violation 
is maximized when Bvz — $13 since J is proportional to the product of Si 2 an d S13. Figure |3]b 
shows R^r for unequal 9\2 and $13: sin 2 2#i 2 = 0.0336 and sin 2 29i 3 = 0.0038, which for 
5 = gives the same LSND result as sin 2 2# 12 = sin 2 2# 12 = 0.015. The CP-violation 
effects for # 12 7^ #13 are not as dramatic as with # 12 = #13, but still may be observable with 
10 20 kt-decays. 

An examination of Table | shows that the relative size of the CP violation in the v e — > 
or v e — > v T channels increases with L/E. However, once A atm is of order unity or larger 
the sin2A atm term in S averages to zero, washing out the CP violation; this does not 
happen until L is much larger than 100 km. Since the flux falls off like 1/L 2 , the statistical 
uncertainty increases roughly like L, and as long as A atm < 1 a wide range of distances have 
comparable sensitivity to CP violation in the v e — > or v e — > v T channels. 

Figure [|a shows the ratio R efl of v e — > event rates (from /i + decays) to v e — * 
event rates (from fi~ decays) versus baseline for oscillation parameters given by Eq. (p"3"D. 
with E^ = 20 GeV. The statistical errors correspond to 2 x 10 21 kt-decays (which could be 
obtained, for example, by three years of running with 10 20 useful muon decays per year, and 
a 10 kt detector having a 67% muon efficiency). A 4 GeV minimum energy cut has been 
made on the detected muon. Although R efl is not as sensitive to CPV effects as R^ T , the 
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increased statistics (due to a larger overall rate resulting from the use of a larger detector for 
muons) make v e — > another attractive channel for CP violation. Figure |]b shows similar 
results for R eT , the ratio of u e —>■ u T event rates (from /i + decays) to v e — > v T event rates 
(from /i~ decays). The statistical errors correspond to 10 20 kt-decays. It is evident from the 
figure that the v e — > v T channel is not as useful for detecting CP violation (primarily because 
of the reduced event rate in tau detection), although this channel is the most sensitive to 
the v e — > v T oscillation amplitude. The combination of parameters IS12S23 — si3C23e l<5 | in the 
v e — > u T amplitude is also present in the sin 2 Alsnd term of Piy^ — > u T ) [see Eq. (0)], along 
with an additional factor involving small mixing angles; hence, it would be more difficult to 
measure |si 2 s 2 3 — Si3C 23 e ll5 | in the — > u T channel. 

Fig. |5] shows, for various CP- violating cases with = 20 GeV, the statistical signifi- 
cance (number of standard deviations) that the ratios R^ T1 Refi, and R eT deviate from their 
expected values for the CP-conserving case. While both the v e — > and — > v T channels 
provide good sensitivity near the optimal L, v e — > is more sensitive for a wider range of 
L, especially for values of 6 that do not give maximal CP violation. An additional potential 
advantage of the v e — > channel is that, in principle, lower energy neutrino factories can 
be used since there is no need to be above the tau-lepton production threshold. However, 
as the energy of muons in a neutrino factory decreases, the sensitivity to CP violation also 
decreases (see Fig. H); this is especially true if a lower bound is imposed on the energy of 
the detected muon. 

In principle, measurements can be made with varying L/E, either by using experiments 
at more than one baseline, or by using a measure of the neutrino energy (for example, 
the total observed event energy) at a fixed baseline. However, since there are only five 
independent quantities in the leading- and subleading-order probabilities in Eqs. (p~0|) — (P~2|) ; 
such measurements still cannot completely determine the three- neutrino parameter set; this 
can only be done at short baselines by a measurement of the subsubleading terms in the 
off-diagonal probabilities, or of the subleading terms in the diagonal — > probability 
[see Eq fllT])], both of which would be very challenging experimentally. Similar conclusions 
apply for additional measurements involving electron detection. 

To summarize, when used in conjunction with long-baseline measurements from K2K, 
MINOS, ICANOE, and OPERA, and with results from MiniBooNE, all six of the three- 
neutrino oscillation parameters can in principle be determined with short-baseline mea- 
surements at a neutrino factory. The short-baseline measurements would determine two 
parameters not determined by other accelerator-based experiments, and would provide a 
consistency check by independently measuring three other parameters also measured by 
other accelerator-based experiments. The v e — > v T channel is most sensitive to the quantity 
IS12S23 — Si3C 2 3e tS \. The v e — > channel provides good sensitivity to CP violation over a 
wide range of L when a large muon detector is used, and the — > v T channel is most useful 
for detecting CP violation near the optimal L for the parameter choice illustrated. 

III. FOUR-NEUTRINO MODELS 
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A. Oscillation formalism 



Four-neutrino models are required to completely describe the solar, atmospheric, and 
LSND data in terms of oscillations, since a third independent mass-squared difference is 
necessary. A fourth neutrino must be sterile, i.e., have negligible interactions, since only 
three neutrinos are measured in Z — > vv decays [[H]]. Following Ref. [p0[| , we label the 
fourth mass eigenvalue itlq. Given the pattern of masses mi, m.2, and from the three 
neutrino case in Sec. II, there is a preferred choice for the scale of mo that can fit all of 
the data, including constraints from accelerator experiments, namely, mo must be nearly 
degenerate with mi, so that there are two pairs of nearly degenerate states separated by 
a mass gap of about 1 eV ||[7|]; see Fig. |l]b. Then 8m\ Q governs the oscillation of solar 
neutrinos, 8m\ 2 governs the oscillation of atmospheric neutrinos, and Sm^i, 8m\ x , 8m\ Q , and 
Sml all contribute to the LSND oscillations. 

Six mixing angles and three (six) phases are needed to parameterize the mixing of four 
Dirac (Majorana) neutrinos; only three of these phases are measurable in neutrino oscil- 
lations. Thus three additional mixing angles, which we label #01, #02, and 6*03, and two 
additional phases are required in extending the three-neutrino phenomenology to the four 
neutrino case. The simplest situation, which occurs in most explicit four-neutrino models, 
is that large mixing occurs only between the nearly degenerate pairs; then the four-neutrino 
mixing matrix can be parametrized as 



U 



( 



Coi 



"Sol 



'01 



Coi 



Coi(.S23 s 02 — C23S03J ■S 01 ( l S23' S o2C23So3j 
-S l(s 2 3-Sl2 - C23S13) +C01O23S12 _ C 23-Sl3) 



V 



'02 



'12 



~~ c 0l( s 23 S 03 + C23S02) — s 0l( S 23 S 03 + c 23 s 02) C 2 3 
+ Soi(s2 3 Si 3 + C23S12) -C01O2V13 + C23S12) 



^03 > 



'13 



'23 



"•§23 c 23 



/ 



(20) 



where Sjk is here defined as sin Oj^e 18 ^ , and the 8jk are the six possible phases for Majorana 
neutrinos. We set 812 = £23 = #02 = without loss of generality, since only three phases are 
measurable in neutrino oscillations. 

In this four-neutrino scenario, the parameters 8m\ 2 -, 8m} ix , 823, Q\2 and #13 have the same 
roles as in the three-neutrino scenario in the previous section; the phase £13 can be identified 
with 8 in the three-neutrino case. The angle #01 describes the mixing of the fourth flavor 
of neutrino with the state it is nearly degenerate with, u e ; together 8m\ and #01 can t & ke 
on the approximate values appropriate to any of the solar neutrino oscillation solutions |32j 
(small angle MSW, large angle MSW, LOW, or vacuum). The remaining mixing angles #02 
and #03 describe mixing of the fourth neutrino with the two states u^, v T in the other nearly- 
degenerate pair. Although models with pure oscillations to sterile neutrinos are disfavored for 



the solar |33] and atmospheric |34j] data, models with mixed oscillations to sterile and active 
neutrinos |7| j35|j36|| , i.e., non-negligible #02 and 9 03 , are presumably also possible. Finally, <5 01 
and 803 are extra phases that may be observable in neutrino oscillations with four neutrinos. 
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In the limit that terms involving the solar mass-squared difference can be ignored, the 
expressions for the oscillation probabilities in Eqs. (10) and (|TT|) remain the same. However, 
the i/p —>■ v T probability becomes 



P{vi* -> v T ) - 4|(si 2 c 2 3 + Si 3 S23e J<5l3 )(si 2 S23 - s 13 c 23 e~ v) 

+ (s 02 C2 3 + s 3S23e i5o3 )(s 02 S23 - s 3C2 3 e~ iSo3 )\ 2 sin 2 A L SND 

+ sin 2 2#2 3 sin 2 A atm + (si 2 si3 sin 5 13 + s 02 s 03 sin S Q3 ) sin 29 23 S . (21) 
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For P{Vfj, — > u T ), the sign of the S-term is reversed. 

B. Discussion 

When combined with measurements from K2K, MINOS, ICANOE, OPERA, and Mini- 
BooNE, the three parameters #12, #13 and 5 13 can in principle be determined by short-baseline 
measurements of the v e — > v T amplitude and the CPV term in the v e — > channel. This is 
similar to the three-neutrino case in Sec. II. However, measurements of — > v r and — > v T 
will only give partial information on 6*02, #03, and 5o 3 . The —>■ survival probability can 
also be measured, which in the limit that the solar mass-squared difference can be ignored 
is 

P(^ -»■ ~ 1 - 4A(1 - A) sin 2 Alsnd - sin 2 26 23 sin 2 A atm , (22) 

where 

A~4 sl 3 (\s Q3 \ 2 + \s 13 \ 2 ) + cl 3 (sl 2 + s 2 12 ) + 2s 23 c 23 Re(so2So 3 + s 12 s 13 ) . (23) 

Both oscillatory terms in Eq. (E^) are second order in small quantities in short-baseline 
experiments. In principle Piv^ — > Vy) could be used as an additional measurement; however, 
unless 6*02 and 6*03 are larger than 12 and #13, the deviations of Piy^ — > v^) from unity are 
very small, and the normalization of the beam flux would have to be known to high precision 
for this to be useful. 

The parameters s\ 2 , s\ 3 and s$ determined by short baseline measurements of the u e — > v T 
and v e — > Ufj, channels give predictions for the — > v T channel at short baselines that can 
differ for the three- and four- neutrino cases [Eqs. (0) and (|2T|), respectively]. If the three- 
neutrino predictions for — > v T are found to substantially disagree with the experimental 
measurements, then the disagreement would provide evidence for the existence of a fourth 
neutrino. The absence of such a disagreement would indicate that either there are only three 
neutrinos or that #02 and #03 are significantly smaller than #12 and 9\ 3 . 

In general, the sensitivity of the measurements of — > v T and — > v T to 6*02, #03 , 
and <5o3 in the four-neutrino case are similar to the sensitivities to d\ 2 , 8\ 3 , and 8 in the 
three-neutrino case. For example, given the parameters in Eq. ( [IB] ) and 5 13 = (i.e., no 
CP violation in the v e —>■ and v e — > v T channels), if S02 = S12 and S03 = S13, then the 
four-neutrino predictions would be given by Figs. ^|a and ||b, where <5 03 takes on the values 
of S in the figures; the corresponding three-neutrino predictions would be given by the 6 = 
curves. Larger values of S02 and S03 could give a much larger CP-violating effect, provided 
that Sq 3 was not small. If both 5q 3 and 5q2 were nonzero, their CP-violating effects could 
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add together either constructively or destructively. Hence, larger CP- violating effects are 
possible in the — > v T channel in the four-neutrino case than with three neutrinos, or 
CP violation may be present in —>■ v T when it is not present in the v e — * or v e — > v T 
channels (or vice versa), unlike the three- neutrino illustrated in Fig. [7]. Even if there 

is no CP violation, effects of the angles 9 02 and 9 03 could be seen in the sin 2 A L snd term 
in Eq. (f2lD , as illustrated in Fig. ||. Here, sensitivities to 6*02 and 6*03 tend to increase with 
decreasing L, due to the increased flux at shorter distances. 

Proof of the existence of a fourth neutrino does not exclude the possibility that there may 
be more than four neutrinos. Strictly speaking, an inconsistency between the measurement 
of — > v T oscillations and the three-neutrino predictions would imply only that there are 
four or more neutrinos. In a model with four or more neutrinos, there are many more mixing 
angles and phases in the neutrino mixing matrix, and it would not be possible to determine 
them all from these oscillation measurements. However, when combined with results at 
other baselines, most of the four-neutrino parameter set could be determined pjf . 

IV. SUMMARY 

In scenarios designed to describe the LSND oscillation results, our results show that 
short-baseline neutrino factory measurements can in principle determine all of the three- 
neutrino parameters provided the results are used together with future results from other 
accelerator-based experiments. The v e — > v T channel is most sensitive to one combination of 
parameters in the CP-conserving terms. The v e — > channel provides good sensitivity to 
CP violation over a wide range of L (20-100 km for = 20 GeV), assuming that a large 
muon detector is used, e.g. 10 kt. The — > v T channel is also sensitive to CP violation 
for a more restricted range of L, and may be used to explore whether more than three 
neutrinos exist. For four or more neutrinos, the CP-violating effect in — > v T may be 
either enhanced or reduced by the additional mixing parameters. If MiniBooNE confirms 
the LSND oscillation results, then it will be important to measure the rates in all three 
appearance modes, as well as to search for CP violation in v e — > and —>■ u T , to obtain 
indirect evidence for the existence of sterile neutrinos. Neutral current measurements would 
complement the charge current studies of this paper in the search for sterile neutrinos. 
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FIG. 1. Ordering and separation of mass eigenvalues in the (a) three-neutrino and 
(b) four-neutrino scenarios in this paper. 



0.9 
0.7 

o 0.5 
" OA 



<v 

i- 

o 

-I-' 
CO 



CO 
CD 

o 
Of 



0.3 

8 
4 

2 
1 



_L (km) 



L (km) 




sin 2 2tf 12 = 0.015 
sin 2 2tf 13 = 0.015 
<5m 2 „ = 0.3 eV 2 



5m z „ = 3.5x1 0" 3 eV : 




-80 -40 40 
<5 (degrees) 



80 



FIG. 2. Ratio of tau production rates for stored fi + to that with stored [i for v e 
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Vfx — ► oscillations, shown versus the CP-violating phase 5 for several values of L. These results 



assume a stored muon energy = 20 GeV and oscillation parameters given by Eq. (13). The 
representative errors shown are statistical, assuming 10 20 kt-decays (after accounting for detector 
efficiency) . 
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FIG. 3. (a) Rut = N(vu — > i> T )/N{v^ — > u T ) versus L, for En = 20 GeV, with oscillation 
parameters given by Eq. (Jl3[) , where 612 = #13. (b) Similar results for sin 2 2^12 = 0.0336 and 
sin 2 2#i3 = 0.0038. Statistical errors correspond to 10 20 kt-decays. 
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FIG. 4. (a) R eu = N{u e -► vJ/Nfa -> and (b) R eT = N(u e -► u T )/N{u e -> P T ). A 4 GeV 
cut has been imposed on the detected muon in (a). Statistical errors correspond to 2x 10 21 kt-decays 
in (a) and 10 20 kt-decays in (b). 
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and 2 x 10 21 kt-decays for muon detection. The other oscillation parameters are given in Eq. (|13|). 
Similar results are obtained for negative values of c5. 
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compared to the CP-conserving case, assuming 2 x 10 21 kt-decays, for = 20 GeV (solid curve), 
10 GeV (dashed), and 5 GeV (dotted). The other oscillation parameters are given in Eq. (^). A 
4 GeV cut has been imposed on the detected muon. 
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FIG. 7. CP-violation effects with four neutrinos in R^ T versus L, for = 20 GeV with 
oscillation parameters given by Eq. (|l3|), O 2 = #03 = 012 = 013, ^03 = (solid curves), 90° 
(dashed), and —90° (dotted) and (a) 5\z = 90° and (b) 8%z = —90°. Statistical errors correspond 
to 10 20 kt-decays. Corresponding results for three neutrinos are given in Fig. |||. 
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